NUB-3263 
YITP-SB-09-10 



O 

o 

(N 



On 



in 
o 



Gluino NLSP, Dark Matter via Gluino 
Coannihilation, and LHC Signatures 

Daniel Feldman 1 , Zuowei Liu 2 , and Pran Nath 1 

1 Department of Physics, Northeastern University, Boston, MA 02115, USA, 
2 C. N. Yang Institute for Theoretical Physics, 



r- 1 . 

o . Stony Brook University, Stony Brook, NY 11794, USA. 

or 

Abstract 

(N 
> 

00 . 

The possibility that the gluino is the next to the lightest supersymmetric particle (NLSP) 
is discussed and it is shown that this situation arises in nonuniversal SUGRA models within 
a significant part of the parameter space compatible with all known experimental bounds. 
It is then shown that the gluino NLSP (GNLSP) models lead to a compressed sfermion 
spectrum with the sleptons often heavier than the squarks at least for the first two gener- 
ations. The relic density here is governed by gluino coannihilation which is responsible for 
a relatively small mass splitting between the gluino and the neutralino masses. Thus the 
GNLSP class of models is very predictive first because the SUSY production cross sections 
at the LHC are dominated by gluino production and second because the gluino production 
itself proceeds dominantly through a single channel which allows for a direct determination 
of the gluino mass and an indirect determination of the neutralino mass due to a linear re- 
lation between these two masses which is highly constrained by coannihilation. A detailed 
analysis of these models shows that the jet production and tagged b-jets from the gluino 
production can be discriminated from the standard model background with appropriate 
cuts. It is found that the GNLSP models can be tested with just 10 fb^ 1 of integrated 
luminosity and may therefore be checked with low luminosity runs in the first data at the 
LHC. Thus if a GNLSP model is realized, the LHC will turn into a gluino factory through a 
profuse production of gluinos with typically only a small fraction < 5% of total SUSY events 
arising from other production modes over the allowed GNLSP model parameter space. 



1 Introduction 



One of the interesting possibilities that arises within the landscape of possible sparticle mass 
hierarchies pQ is that the gluino (g) is the next to the lightest supersymmetric particle (NLSP) 
where neutralino dark matter produces the correct relic abundance of such matter consistent 
with the WMAP observations [2]. In fact, an analysis in the context of nonuniversal supergravity 
models (NUSUGRA) reveals that the gluino NLSP model (GNLSP) arises in a significant part 
of the parameter space (3j Hj. Amongst the various possible ways that the first four lightest 
sparticles may stack up in their mass hierarchy, one finds three such hierarchical mass patterns 
where the gluino is the NLSP [HI Hj which have been classified as models (NUSP13, NUSP14, 
NUSP15) ETas given in Table®. We will often refer to this subclass of NUSUGRA as the GNLSP 
class of models. Although progress has been made on the parameter space of sparticle masses 
with coannhilating gluinos [SI El El IU Cj , collider and dark matter detection implications of the 
GNLSP models have yet to be explored in any great detail. Thus in this work we give a dedicated 
analysis of such a model. Since the gluino is a strongly interacting particle, an NLSP gluino will 
change drastically the typical sparticle analyses. The phenomenology of GNLSP models is very 
different from that of a model where the gluino is the LSP [8] which we do not discuss in this 
paper. We note also that while relatively light gluinos have been studied in detail in reduced 
SU(3) gaugino mass models [9], the GNLSP situation, which we cover here, was not explored. 



Table 1: Hierarchical sparticle mass patterns for the four lightest sparticles, where x = Xi is 
the LSP neutralino, and where the gluino is the NLSP that arises in the NUSUGRA models. 
The labeling of the mass patterns is as given in [HI H] . 



The outline of the rest of the paper is as follows: In Sec© we discuss the origin of nonuni- 
versalities in the gaugino masses in 577(5), 5*0(10) and E s GUT models. Thus this means that 
the ratio of the U(l), SU(2) and SU(3) gaugino masses at the grand unification scale Mq are not 
in the ratio 1:1:1. Here we point out that while no F term breaking with a single irreducible 
representation can generate a GNLSP model, it is possible to do so with a mixture of two (or 
more) such breakings. Specifically we consider a linear combination of breaking with a singlet 
and a non-singlet F term and show that several models exist which lead to a GNLSP model. We 
also show that there exists a subclass of models which superficially look different but are in fact 
isomorphic. In Sec.([H]) we discuss the techniques for the computation of the sparticle spectrum at 

1 There is another sparticle mass pattern NUSP10 [I] with the mass hierarchy x° < t% < g < xf, where t\ is 
the NLSP, but the g lies close to the t\ mass. 



NUSP 



Mass Pattern 



NUSP13 
NUSP14 
NUSP15 



X° < 9 < Xf <X% 

x° < g < h < xf 
x° <g <A~H 
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the weak scale and also discuss the experimental constraints that are imposed on the spectrum. 
In Sec. ((!]) we give a discussion of how the relic density consistent with WMAP data is satisfied 
under the assumption that dark matter is entirely constituted of cold dark matter in the form of 
R parity odd LSP neutralinos. Here it is shown that there are two main mechanisms by which 
this can come about. The first mechanism is by coannihilation with gluinos where the dominant 
processes which participate in the coannihilation are x°x° ~~ * f f i X°9 99 ~~ * 99i QQ- The 

second mechanism is the one where the LSP has a significant higgsino component and here the 
relic density constraint is satisfied in a similar fashion as in the usual higgsino dominated LSP 
model. However, this is rather rare, and when it occurs, it is often with a small amount of gluino 
coannihilation. In Sec. (jSJ) we delineate the allowed parameter space for the GNLSP models and 
show that there is a significant region of the parameter space where such models manifest. In 
Sec. (El) we show that the GNLSP models lead to a compressed sfermion spectrum for the first two 
generations. Specifically the sleptons and the squarks of the first two generations are essentially 
degenerate with the sleptons sometimes being heavier than the squarks. In Sec. (CO) we give an 
analysis of the signatures for GNLSP models. Here we discuss three sets of post trigger level cuts 
labeled CI, C2 and C3 which are designed to reduce the background and enhance the signal to 
background ratio. It is found that the dominant signatures are jets and missing energy and with 
properly chosen post trigger level cuts they stand out above the background. It is further found 
that a discovery of a GNLSP model can come about with an integrated luminosity of 10 fb _1 at 
the LHC with a gluino mass of up to 800 GeV. More generally, the models discussed here can 
be put to test with the first data from the LHC. 

In Sec. OH]) we discuss the direct detection of dark matter in GNLSP models. It is found that 
the CDMS-08 data already constrains the parameter space of GNLSP models, although only 
rather mildly at the level of &si{x°p) — 10 _44 cm 2 . Further, the future data from CDMS and 
LUX will either detect dark matter predicted in this model or constrain the parameter space of 
the model. It is also noted, however, that a part of the parameter space of the model leads to a 
rather small spin independent neutralino-proton cross sections, i.e., &si{x°p) < 10 _46 cm 2 , which 
lies outside the reach of the current direct detection experiments and similar experiments in the 
foreseeable future. Interestingly, much of this parameter space will be accessible at the LHC since 
the gluinos can be produced and detected via their jet and missing energy signatures as discussed 
in Sec.(J7]). In this sense, the LHC and the direct detection experiments are complementary. In 
Sec. ([9]) we discuss the benchmarks for the three GNLSP model sets A, B and C. 

Conclusions are given in Sec. (ITOl) . In Appendix A we give sum rules on the gaugino masses 
that hold for the various cases of nonuniversalities that appear in Table ([2]). These sum rules also 
hold when one includes a singlet breaking along with breaking with a non-singlet. In Appendix 
B we give an analysis at one loop which explains the compression of the sfermion spectrum for 
the first two generations. Benchmarks discussed in Sec. (Q are given in Appendix C. 
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2 Gaugino mass nonuniversalities in GUT models, the 
gluino NLSP and scaling 



There is considerable literature on nonuniversalities of soft breaking and their applications |10] 
within the framework of supergravity grand unification [TTI IT2"] . Our focus will be on the gluino 
phenomenology that results from the gluino being an NLSP [For recent analyses related to gluino 
phenomenology in various contexts see [H]]. Specifically, our focus here will be on nonuniver- 
salities in the gaugino mass sector arising from F type breaking in SU(5), SO (10), and E 6 GUT 
groups which have been discussed over the years [TH [15] and a more comprehensive analysis has 
been given recently [16]. Results of this analysis are summarized in Table© H. In the table, 
ratios of gaugino masses that arise when the GUT symmetry is broken by an F term, which is an 
irreducible representation of the gauge group SU(5), 5*0(10), and E 6 , and enters in the decom- 
position of the symmetric product of two adjoint representations corresponding to the relevant 
group. Table© identifies the group and the irreducible representation and the corresponding 
ratio of the gaugino masses. For SO (10) and E s several gaugino mass ratios are listed for a given 
irreducible representation. These correspond to different patterns by which the GUT symmetry 
breaks to lower rank groups. Further details can be found in [16] . None of the models listed in 



Group 


Rep. 


Label 


Mi : M 2 : M 3 


Group 


Rep. 


Label 


Af i : M 2 : M 3 


SU(5) 


1 




1:1:1 




650 


(12) 


-1:1:1 




24 


(1) 


-l/2:-3/2 : 1 






(13) 


-1:1:0 




75 


(2) 


-5:3:1 






(14) 


1/10 : -3/2 : 1 




200 


(3) 


10 : 2 : 1 






(15) 


-13/5 : 1 : 1 


50(10) 


210 


(4) 


-3/5 : 1 : 






(16) 


1/5:1:0 






(5) 


-4/5 : : 1 






(17) 


41/15 : 1 : 1 






(6) 


1:0:0 




2430 


(18) 


-11/5 : 1 : 




770 


(7) 


19/10 : 5/2 : 1 






(19) 


1 : 35/9 : 1 






(8) 


32/5 : : 






(20) 


12/5 : : 1 




650 


(9) 


-1/5 : 1 : 






(21) 


0:0:1 






(10) 


-1/5 : -1 : 1 






(22) 


33/5 : 1 : 1 






(11) 


3:1:1 






(23) 


9/5:1:0 



Table 2: Exhibition of the gaugino mass ratios at the GUT scale for various groups and represen- 
tations in SU(5), 5O(10), and E 6 modelsfTH]. The mass ratios are listed in a hierarchical manner, 
i.e., they are listed in the order of the smallest rank group and lowest dimensional representation 
in which they first appear and are labeled from (l)-(23). Thus a specific ratio may be repeated 
several times as one goes up the chain. 

Table© can give rise to a gluino as the NLSP with F type breaking with a single irreducible 
2 In this analysis we do not consider the flipped models and the ratios listed in Table ((2]) exclude such models. 
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represent at ionE). However, we will show that a combination of GUT symmetry breaking in the 
gaugino mass term sector with two irreducible representations does allow for a gluino as the 
NLSP for a subset of models listed in Table([2]). Specifically we will consider a linear combination 
of a singlet and a non-singlet F term. In this case an interesting phenomenon arises in that the 
models with the same value r = (M 2 — Mi)/ (Ms — Mi) can be made isomorphic under redefini- 
tions and scalings in the gaugino sector. Thus suppose we write the gaugino masses for models 
of the above type with a singlet and a non-singlet F breaking so that 

Aff = (1 + aiai)m 1/2 , M 2 W = (1 + b^m^, M 3 (i) = (1 + c^)m 1/2 , (1) 

where the first term within each of the parentheses on the right hand side in Eq. (prj arises from 
the singlet contribution, and the second term within each of the parentheses is the contribution 
from the non-singlet. Here i defines a specific model and 6j, c, are the fractions given in 
Tabled with a, being an arbitrary parameter. Next we note that two models % and j defined 
by Eq. ([T|) can be made isomorphic if they have the same value of r in the sense that 

MP = \ ij MP;a = 1,2,3, (2) 

when atj is related to a, in the following way 

a^ l (h - di) = a^ibj - aj) + ajfy - (3) 

This means that under the constraint of Eq.([3]), a rescaling of m 1 / 2 of model j can make it 
isomorphic to model %. Thus in essence, models with the same value of r would in fact be 
equivalent when taken in a linear combination of breakings including singlets. Using Eq. ([1]) and 
Table([2]) one finds that there are several possibilities for which the GNLSP class of models can 
arise. We limit ourselves to the following cases: 



1. Model GNLSPa (ISO-I) : This class of models arise where r takes the common value —2/3 
as exhibited below 

Mi : M 2 : M 3 
-1/2 : -3/2 : 1 
19/10 : 5/2 : 1 
-1/5 : -1 : 1 

2. Model GNLSPb : This is an E 6 model with F type breaking with 2430 plet such that [15] 
E 6 -> 377(6)" x ,S77(2) L (2430 -> (189, 1)) which gives M x : M 2 : M 3 = : : 1. This model 
can generate a gluino as the NLSP upon the addition of breaking with a sing letfl 

3 This also holds for flipped models, i.e., F type breaking with a single irreducible representation cannot give 
rise to a gluino as the NLSP. 

4 We note that there is another Eq model with F type breaking with 2430 plet such that Eq — > SU(6)" x 
5f/(2) L (2430 ->■ (405, 1))[S] which gives Mi : M 2 : M 3 = ^ : : 1 (r = 12/7). This model can also generate a 
gluino as the NLSP upon addition of a singlet and there is a relative sign flip between Mi and M 2 in this case. 
However, the model gives a light Higgs mass in the parameter space investigated which falls below the current 
limits and thus we do not consider this model further. 
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3. Model GNLSPc : Here r is free and thus defining r = 62/63 the gaugino masses at the 
GUT scale may be parametrized as 

M 1 =m 1/2 , M 2 = (1 + 6 2 )m 1/2 , M 3 = (1 + 6 3 )m 1/2 , (5) 

and 6 2 and 63 can be varied independently. Model GNLSPc contains models GNLSP^ and 
GNLSPb as subcases^]. 

Aside from the model discussed in footnote 4, models GNLSPa , GNLSPb , and GNLSPc are the 
only models which lead to a GNLSP through breaking with a singlet and a nonsinglet. This can 
be seen easily by using the semi analytic analysis given in Appendix B. For all the three models 
a GNLSP requires 63 to lie in the range (—0.9,-0.8). Some benchmarks for Models A,B and 
C are given in Tables (loTlSTTUj) and a display of their partial sparticle spectrum and some other 
properties of these models are exhibited in Tables (!Tll9l[TTj) . We also note that from the analysis 
of [TB] one can discern another set of models which have the same common value of r. Thus the 
models with the gaugino mass ratios Mi : M 2 : M 3 = (i) — ~ : 3 : 1; (ii)| : 2 : 1; (hi) — | : 1 : 0;(iv) 
I : — I : 1; (v) ^ : | : l;(vi) | : : 1 have the common value r = 8/3. One may call this ISO-II 
because when combined with a singlet F type breaking these models too would be isomorphic so 
that the six different models are effectively one model as far as the gaugino sector is concerned. 
However, this model class does not lead to a GNLSP which is the main focus of this paper. In 
the following we discuss the GNLSP models in further detail including the satisfaction of the 
relic density, the production cross section of the gluinos, the signatures for their identification at 
the LHC, and the direct detection of dark matter in the GNLSP class of models. 

3 Experimental constraints 

Our general procedure is similar to that discussed in [4] which we briefly describe below. In 
the analysis one specifies boundary conditions of the model at the grand unification (GUT) 
scale which we take to be Mq ~2x 10 16 GeV. Specifically we take the sfermion masses at the 
GUT scale to be universal, but assume that the gaugino masses are in general nonuniversal with 
nonuniversalities given by 6 2 and ($3. One then uses renormalization group equations (RGEs) to 
compute the sparticle mass matrices and their eigenvalues at the electroweak scale. The code 
used in these RGE evolutions and computations of the sparticle spectrum is SuSpect2.41 [IT], and 
similar results are obtained with SoftSUSY[18] and SPheno [19]. Further, one imposes the lower 
limit constraints on the sparticle masses from the LEP and from the Tevatron data as well as 
constraints from the WMAP on the relic density. The analysis of the relic density is first done at 
the perturbative level with MicrOMEGAs [20], which relies on CalcHEP [21]. Non-perturbative 
effects on the relic density are also discussed. 

5 We remark that in [H [4] , where the gluino NLSP in SUGRA models was previously observed, the notation 
^5 j $6 was used. 



5 



Below we give a list of the relevant constraints from collider and astrophysical data which 
have been included in the analysis: 

(i) The 5- year WMAP data constrains the relic density of dark matter in the universe so that 
Odm^ 2 = 0.1131 ± 0.0034 [22]. We take a 6a corridor around the central value to constrain the 
relic abundance of neutralinos. The larger band is taken due to the sensitivity of the relic density 
computation in particular regions of the parameter space. A large class of our models fall well 
within a 2a bound. 

(ii) The FCNC process b — > s'j receives a significant contribution from the SUSY processes |23j. 
The Heavy Flavor Averaging Group (HFAG) [24] along with the BABAR, Belle and CLEO give 
experimental results: Br(B -> X s ^) = (352 ± 23 ± 9) x 10~ 6 . A new estimate of standard 
model contributions at 0(a 2 s ) gives [25] Brib — > sj) = (3.15 ± 0.23) x 10~ 4 . We utilize both 
experimental and theoretical progress in the evaluation of this observable and take a 3a corridor 
around the experimental value, 2.77 x 10 -4 < Br(b — > sj) < 4.27 x 10 -4 , to constrain the 
theoretical prediction including both SM and SUSY contributions. 

(iii) Another important constraint from B-physics is the rare decay process B s — > fi + fi~ 
which can become significant for large tan/3[26]. The most stringent 95% (90%) CL. limits are 
achieved by CDF[27] Br(B s -> /i + ^") < 5.8 x 10" 8 (4.7 x 10" 8 ). We take a conservative limit 
Br(B s -> fi+fi-) < 10~ 7 . 

(iv) For the constraints from the anomalous magnetic moment of the muon, we use a con- 
servative bound —11.4 x 10~ 10 < 5(g M — 2) < 9.4 x 10~ 9 as in [28] where 8{g fJi — 2) is the new 
physics contribution to (g M — 2) beyond the standard model. 

(v) Additionally, we also impose various mass limits as follows: m~± > 104.5 GeV [29] for the 
lighter chargino, m-j- > 101.5 GeV for the lighter stop, and > 98.8 GeV for the lighter stau. 
For the lightest CP even Higgs boson mass in MSSM we take the constraint to be > 100 GeV 
(90% of the models that pass all constraints have > 110 GeV). One may compare these with 
the standard model like Higgs boson mass limit which is ~ 114.4 GeV [30]. For the gluino mass, 
recent Tevatron experiments give m g > 308 GeV (D-Zero)[3T] and m g > 280 GeV (CDF) [32] . 
The limits given by [3TL [32] are valid within the framework of the minimal supergravity models 
and may be modified in nonuniversal SUGRA models. Hence the total SUSY production cross 
section constrained by the Tevatron analyses will typically be a larger total cross section than 
that which arises in the GNLSP models. Further, as we will show shortly, the mass splitting 
between the NLSP gluino and LSP neutralino must be relatively small in order to satisfy relic 
density constraints. Thus the relatively small mass splittings between the LSP and GNSLP can 
lead to softer decay products and an overall lower multiplicity of final state events relative to 
models for which the mass splitting is significantly larger. Therefore in this analysis, we take 
a conservative lower bound, namely m g > 220 GeV. Our choice of this lower bound is taken 
as to not eliminate a part of the parameter space which may otherwise be allowed pending a 
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full analysis of the Tevatron data using nonuniversalities (see also [33] for a related discussion 
regarding a lower bound on the mass of the gluino). 



4 Relic density via gluino coannihilation 

It is interesting to ask how the relic density constraints are satisfied in the class of models with the 
gluino as the NLSP as these constraints have important implications for collider phenomenology 
(for recent works connecting sparticle phenomenology at colliders and dark matter see [3U 
I3T)1 I3"7] 135]). As an illustration we consider the model 

GNLSPco : (m ,m 1/2 , A , tan p,S 2 ,S 3 ) = (1450,730,2700,40,0.332,-0.839), (6) 

where all masses are in GeV and sign(fi) is taken to be positive. We will take the top mass at 
170.9 GeV throughout this work, though the analysis here does not show great sensitivity to the 
top mass. This model gives (m^,mj) = (305.1,348.6) GeV. For GNLSPco the channels which 
contribute to l/(Qh 2 )^o more than 1% are as follows : gg — > gg(47%), gg — > uu(8%), gg — > 
cc(8%), gg - dd{8%), g~g - ss(8%), ~gg - 66(6%), gg - tt(4%), x°x° - 66(6%), x°~g - tt(2%), 
X°X° -> tt(2%), x°X° -> r+T - (l%). The relic density is (fi/i 2 )^o = 0.108 at the perturbative 
level, and the model has eigen decomposition x° — 0.9866 — 0.016tZ> + 0.146/ii — 0.092/i2 where 
b,w, are the bino and wino components and hi,fi2 are the higgsino components, and thus the 
model has a substantial higgsino component. This model belongs to the pattern classified as 
NUSP13 in[3lig. 



From the above it is clear that the gluino processes dominate the WIMP annihilation at 
the freezeout temperature in the early universe Further the LSP mass and the NLSP mass 
are close with a mass difference A^o = (m§ — m^o)/m^o m 0.14. An examination of the mass 
splittings and the associated annihilation processes point to a strong coannihilation occurring 
in the model of Eq.([6"l). Thus consider the annihilation processes XiXj g°i n g into the standard 
model particles. Here the effects of coannihilation are controlled by the Boltzmann suppression 
fact or [39] 

where gi are the degrees of freedom of Xi, x — mi/T and Aj = (m; — mi) /mi, with mi defined 
as the LSP mass. The processes which dominate the WIMP annihilation in the early universe 
are 

X V - F, x°~9 - F', ~gg - F" , (8) 

where F, F', F" constitute the pairs of standard model states. The relic density is controlled by 
the integral 

POO 

J X} = / x~ 2 (<j cS v}dx , (9) 

J X f 
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where v is the relative velocity of annihilating supersymmetric particles, (cr c sv) is the thermally 
averaged cross section times the relative velocity and Xf is the freezeout temperature. The cr e g 
that enters the relic density can be written approximately as follows 



* o&fa 7 2 + 2 7 ^ + ^2c \ (10) 
V a gg a gg J 

where 7 = ig/lx an ^ wnere 7« are defined by Eq.(j7J) and where [SJ 

'($9 -+99) = ^{logi±|[21-6/3 2 -3/3 4 ]-33/5 + 17/3 3 }, 



Here /3 = yl — 4m?/s, and the quark mass enters Eq.l JTTl) through /3 = yl — 4m^/s. One 
interesting phenomenon concerns the following: we know that the cross section for the annihi- 
lating gluinos falls with the gluino mass. On the other hand a c s that enters the relic density 
analysis must be nearly constant for a wide range of gluino masses so that the relic density be 
satisfied. This can happen due to the presence of the coannihilation factor 7 2 which multiplies 
a^g in Eq. ffTOl) . This is easily seen by noticing that A^o = {m g — m x o)/m x o has a dependence 
on the gluino mass of the form 

AjL„ - C log(m~ g /ml) , (12) 

where A9~ = (m~ — m°~ )/m°~ , and m~ is a reference gluino mass and m°~ is the corresponding 
reference neutralino mass, and C > 0. What one finds is that the difference A^o decreases when 
gluino mass increases which enhances 7 and compensates for the falling cross section a gg . The 
above phenomenon sustains an essentially constant J Xf as the gluino mass varies allowing for a 
satisfaction of the relic density over a wide range of gluino masses. What the analysis implies is 
that rrig/m^o tends to unity as the gluino mass increases. A numerical analysis bears this out. It 
is, however, interesting to note, that we also find few cases where the GNLSP emerges without 
significant coannihilation which occurs when the LSP has a significant higgsino component which 
allows for the satisfaction of the relic density constraint in a manner quite similar to what 
happens on the Hyperbolic Branch of REWSB[40j. Indeed, there are cases where the neutralino 
annihilations are seen to dominate the annihilation cross sections via x°X° ~^ {bb,T + r~) with 
only a small contribution to the satisfaction of the relic density constraints arising from LSP- 
GNLSP coannihilation, and in some cases coannihilation enters only at the single percent level (an 
example is GNLSPci given in Appendix C). Another interesting example is model GNLSPai (also 
given in Appendix C) which proceeds with annihilation contributions to relic density calculation 
dominantly via (62%)x°X° ~* tt, and (15%) x°X° ~^ W + W~, and only a small fraction (3%) for 
gg — > gg and the remainder coming from neutral diboson final states. 

We discuss now possible nonperturbative corrections to the annihilation cross section. As 
shown in Refs [SI [S] nonperturbative effects on the annihilation cross section can be relevant near 



3 A similar relationship in a graphical form appears in the analysis of [5] for a bino LSP. 
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threshold where multiple gluon exchange, for example, can give rise to the so-called Sommerfeld 
enhancement factor £. These effects may be approximated as [8] 





Cjira s 1 


1 — exp | 


- p w 





(13) 



where Cj =g = 1/2 (Cj =q = 3/2) for gg —> gg (gg —> qq) respectively, and we note that Sj enters 
bilinear ly in Eq. liTTj) . 



Bound states can form as well if the gluino is a stable LSP. We do not consider the latter 
situation. However, as already discussed, for the GNLSP at the perturbative level, the dominant 
contribution in the annihilation cross section for most models arises from the gluino-gluino an- 
nihilation modes, and this occurs for A^o < 0.2. Since micrOMEGAs performs the relic density 
analysis using only perturbative cross section, we have carried out an independent analysis of the 
relic density to include the effects of the Sommerfeld enhancement. Our analysis gives results 
which are in agreement with the analysis of third reference of [5] . We note here that an increase 
in Ajj^o in the range of (2 — 3)% is needed when the Sommerfeld enhancement of cross section is 
taken into account. Equivalently the effect of the Sommerfeld enhancement can be recast as a 
shift in the gluino mass for a fixed LSP mass in order to have the same relic density as for the 
perturbative case. Specifically, an upward shift of the gluino mass by a few GeV is needed. Thus 
for example, for the model GNLSPco discussed above, the relic density constraint is satisfied 
with the inclusion of non-perturbative effects with a 3 GeV upward shift of the SU (3) gaugino 
mass at the GUT scale leading to an increase in A^o = 14% — > 16%. More generally, we find 
that numerically for the non-perturbative case, for fixed LSP mass, m g needs to be increased by 
(3 to 6) GeV to achieve the same relic density as for the perturbative case. The above holds for 
gluino masses in the range up to about 1 TeV. In Appendix C we give benchmarks including the 
effects of the Sommerfeld enhancement. 



5 Consistent parameter space of the gluino NLSP 



Based on the initial discovery of the existence of a viable parameter space where the gluino is 
the NLSP in nonuniversal SUGRA models[Hl H], we perform in this work a dedicated search 
for delineating the parameter space of GNLSP models consistent with the radiative electroweak 
symmetry breaking constraints and with the experimental constraints from colliders and from the 
relic density. In our analysis the input parameters m , myz, A , tan /3, <5 2 , #3 assume the following 
bounds: m < 4 TeV, myz < 2 TeV, |Ao/mo| < 3, and tan (3 G (1, 60), while 82 and #3 are chosen 
in a manner appropriate for models A, B and C defined in Sec.(j2])- Thus for model GNLSPa , 
82 is determined by the constraint r = —2/3, for model GNLSPb , 82 — and for the model 
GNLSPc , 82 is assumed to lie in the range 82 G (—0.9, 1), while £3 typically lies in the range 
£3 G (—0.9, —0.8). Within the ranges assumed above, we find a significant region of the parameter 
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Figure 1: (Color online) An analysis of the consistent parameter space in GNLSP models. Red 
(dark) model points are for Model GNLSPc , while the yellow (light) and the green model points 
are for Model GNLSPa and Model GNLSPb as discussed in the text. A e d is the splitting of the 
slepton and squark masses in the first two generations as also discussed in the text. 



space where each model is realized. Of course, the parameter space for model GNLSPc is larger 
than that for the model GNLSPa or for the model GNLSPb , but the parameter space for 
models GNLSPa and GNLSPb are also quite significant as shown in Fig.([T|). Some interesting 
observations can be made from the analysis of Fig.([T]). Thus the top left panel of Fig.([T|) shows 
that typically rriQ > mi/2 for this class of models while the top right panel shows that the region 
Ao/tuq = is very thinly populated which is in sharp contrast to the mSUGRA case where the 
Ao/tuq region is heavily populated. The lower right panel of Fig.([T]) displays the allowed model 
points in the $3 vs 62 plane which shows that GNLSP models constrain the nonuniversality $3 to 
lie in a very narrow range (—0.9, — 0.8) [HI [4] while 62 is widely dispersed for model GNLSPc but 
restricted for models GNLSPa and GNLSPb since 62/63 = —2/3 for model GNLSPa and #2 = 
for model GNLSPb • These points are indicated in yellow (model GNLSPa) and green (model 
GNLSPb). 
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Figure 2: (Color online) Left panel: An exhibition of the scaling between the light chargino mass 
and the LSP mass for the GNLSP models GNLSP A , GNLSP B and GNLSP C vs the mSUGRA 
model. The figure shows that GNLSPa produces the ratio m^±/m^o ~ 3 which differentiates it 
from the bino branch of mSUGRA and for GNLSPb- Right panel: An exhibition of scaling in 
rrig/m^o. All GNLSP models are well separated from mSUGRA in this figure. 



6 Compression of the sparticle mass spectrum in GNLSP 
models 

In models with universal boundary conditions at the GUT scale for the gaugino masses, the 
gluino mass will be typically a factor of 5-6 larger than the lightest neutralino mass. A large 
gluino mass tends to contribute a significant portion to the squark masses in the RG running. 
Thus in the mSUGRA modelpj] there is typically a significant splitting of the slepton and squark 
masses for regions of the parameter space where mo and mi/2 are comparable. This is generally 
not the case in the model under consideration where gluino is the NLSP. Here the gluino mass 
will be typically much lighter relative to the squark masses, and thus splittings between sleptons 
and squarks will be less pronounced. Specifically for the first two generations one will find a 
rather compressed spectrum. Table ([3]) exhibits the high degree of degeneracy of the squarks and 
of the sleptons in the first two generations in models with the gluino as the NLSP relative to, 
e.g., the mSUGRA model [For a general discussion of sum rules see [H]]. In the examples shown 
one finds that while for the mSUGRA m spdll case the splitting between the sum of the down 
type quarks and the charged sleptons in the first generation can be as much as 35%, while for the 
GNLSP case it is only about 1%. Further, while for the mSUGRA case the first and the second 
generation squarks are invariably heavier than their corresponding slepton counterparts, for the 

7 The mSUGRA pattern mSP3 has the following hierarchy for the first four sparticles: \° < xt < X2 < ^l- The 
largeness of the sfermion masses indicates that the electroweak symmetry breaking is realized on the Hyperbolic 
Branch[40] (for recent works on the HB see [42], [37]) ■ 
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GNLSP case one finds that one can often get an inversion, i.e., the model gives rise to sleptons 
heavier than their squark counterparts within the first and second generations. Specifically, 
defining 

n\ (ma -\- Tns ) — (m?, ■+ Trip n .) 

a W _ o V gij d 2 iJ \ e i» e ^> ■ _ i 9 (-1 a\ 

[ m d u + m dJ + ( m eu + mej 

where i is the generation index, for the mSUGRA case one finds that A^J are positive and 
typically a significant fraction. However, for the GNLSP case one has 

|A2|«1, (15) 

and often |A^| lie in the range much smaller than 1%. Thus the validity of Eq. (fT5]) implies a 
high degree of degeneracy of the squark and slepton masses for the first two generations, and the 
observation of such a degeneracy will provide a strong corroborating evidence along with collider 
signals for testing the validity of the GNLSP models. Of course, a test of Eq. ffT5]) would require 
determination of the squark and slepton masses with a certain degree of accuracy. 



More generally the lower left hand panel of Fig.([T]) exhibits A e ^ as function of Ac = 
where the gluino NLSP and neutralino LSP coannihilate to produce the consistent relic density 
observations of WMAP. The analysis of this panel exhibits more generally the results of Table© 
in that one finds that in all these models A e ^ is relatively small and often negative. We thus arrive 
at the important general conclusion that in the model where the gluino is the NLSP one gets a 
compressed sfermion spectrunj§ for the first two generations relative to the squarks, with mass 
differences between squarks and their slepton counterparts which are typically order a few percent 
and often less over a wide range of the parameters space. In Fig. ([2]) the left panel shows that the 
model GNLSPa can be discriminated from the bino branch of mSUGRA, while the right panel 
of Fig. ([2]) shows that all GNLSP models can be discriminated from mSUGRA. Additionally 
a comparison of the right and the left panels also allows a discrimination of GNLSPa from 
GNLSPb • 



Model 


Pattern 


77lei + rrie 2 


m d, + m ~do 


a£ 


mSUGRA 


mSP3 


5377 


7652 


35% 


NUSUGRA SU(5) 


NUSP13 


7386 


7373 


-0.1% 


NUSUGRA SO(10) 


NUSP13 


7369 


7300 


-0.1% 



Table 3: Exhibition of the mass compression for sleptons and squarks in the first two generations 
in two typical GNLSP SUGRA models with nonuniversalities and a comparison with a mSP3 
model point in the mSUGRA case. Mass splittings between the sleptons and the squarks are 
seen to be much smaller for the GNSLP models compared to the mSUGRA case. 

8 The compressed sfermion spectrum discussed here is very different from the one discussed in[43j. 
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7 Signature analysis at the LHC 



Gluino Mass and Gluino Production Cross Sections: The production cross sections for 
gluinos were studied early on [Hj and the NLO evaluations have also been given [35]. A par- 
ticularly interesting situation is the one which is discussed in the preceding sections where the 
gluino is the NLSP, as this possibility leads to a rather predictive model. Thus one finds that 
in the GNLSP case the gluino production cross section dominates all other SUSY processes and 
further the production is controlled to a large degree by a single process which is gg —> gg, i.e., 
(jpp(SUSY) « a pp (gg — > gg) where cr(SUSY) is the LHC production cross section including all 



2^2 SUSY production modesjlS]. A numerical analysis of the above is shown in Fig.Q for 
the GNLSPs. A consequence of the dominance of the processes gg — > gg over all others has 
interesting implications. Specifically it opens up the interesting possibility that a rather precise 
determination of the gluino mass can be made from a measurement of the production cross sec- 
tion of all SUSY processes. Further since the neutralino is linearly related to the gluino through 
Ag^o the LHC data in this case could allow us to determine the neutralino mass with a fair 
degree of accuracy should the gluino mass be reconstructed. 



LHC as a Gluino Factory 
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Figure 3: (Color online) Left panel: Display of pp cross sections including the individual produc- 
tion modes arising from subprocesses gg — > gg and qq — > gg and the total SUSY cross sections 
plotted as a function of the gluino mass showing the dominance of the gg production process. 
The analysis of the figures shows that in the GNLSP case the LHC will turn into a gluino factory. 
Right panel: Total number of SUSY events passing the Level 1 (LI) trigger cuts and post trigger 
level cuts CI as defined in the text. A majority of the events pass the triggers as can be seen 
by comparing with the left and the right panels of Fig. ([3]) and taking into account that the right 
panel of this figure is given for 10 fb _1 of luminosity. 



It is worthwhile to pause and comment on the sensitivity of the relic density to the codes that 
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are often used for its computation. We exhibit this sensitivity in Table (141) where a small dialing 
of parameters has been done to keep the relic density in the corridor allowed by WMAP. The 
analysis shows that the sensitivity to the codes is rather small and the GNLSP model is robust in 
that it appears in all the three codes used in Table(j4]). In Table([5]) we compare the leading order 
(LO) predictions of PYTHIA and PROSPINO in the GNLSP model for the same parameter 
point given in Table(j4]), and show the next to leading order (NLO) prediction using PROSPINO. 
We also compare the relevant branching ratios. One observes that the squarks decay back into 





SuSpect(2.41)jT7] 


SOFTSUSY(3.0.2)[18J 


SPheno(2.2.3)[19j 


x° 


336.3 


334.5 


334.5 


9 


382.7 


379.8 


381.0 


xf 


424.1 


422.4 


422.9 


h 


451.4 


464.7 


447.4 




0.115 


0.105 


0.117 


{62,63) 


(-0.340, -0.835) 


(-0.340, -0.824) 


(-0.340, -0.823) 



Table 4: A comparison of the sparticle spectra and of the relic density for a GNLSP model with 
(mo, mi/2, A , tan (3) = (1387, 792, 3026, 27) (all masses in GeV), with // > for m t = 170.9 GeV 
and 62 = —0.340 with mb(rrib) and a s (Mz) taken with the default values of MicrOMEGAs (MO). 
For these models Xi < X2 (NUSP13). We use MO 2.2CPC for the first two cases and MO 2.07 
for SPheno. This particular comparison is made at the perturbative level. 



PYTHIA[49J a LO 



99 -> 99 = 70 pb 
99 = 6.2 pb 
Qj9 QjL9 = 1-1 pb 
99 -» hh = 1.5 pb 
else <C 1 pb 



m — 

SUSY 



80 pb 



PROSPINO [50] a LO 


-Knlo 


C"NLO 


~g~9 = 84.3 pb 


1.72 


145 pb 


qg = 3.12 pb 


1.60 


5.0 pb 


tJi = 0.80 pb 


1.55 


1.24 pb 


else <C 1 pb 






a^ SY = 88.5 pb 




a^ Y = 151.6 pb 



Decay 


BR PYTHIA 09] 


BR SUSY-HIT [47J 


9 -> (bbx°, uux°, ddx°) 
9 -> X°9 

Ql -> g(u,d) L 

h -> x + b 


(20,61,19)% 

82% 
100% 


(20,61,19)% 
0.03% 
86% 
100% 



Table 5: A specific exhibition of the dominance of the process gg 
and NLO for the GNLSP model point given in TablefllJ). 



gg in pp collisions at LO 



a gluino with a large branching ratio, 
at the level of 96%, i.e., a 



NLO 



The NLO calculation retains the dominance of the o"^ LO 



0.96 a*™ . 



We note that the model in Table© has a very 
small branching fraction g — > x°9 (calculated with SUSY-HIT j4T] and a very similar suppression 
is seen with ISAJET gg] (<0.1%)). This is to be contrasted with the GNLSP Co model (see 
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Figure 4: (Color online) Left panel: A display of the average squark mass vs the gluino mass 
in GNLSPc and a comparison with mSUGRA. Right panel: A display of the lightest squark 
mass vs the gluino mass in GNLSPc and a comparison with mSUGRA. The gluino mass is 
relatively light in the GNLSP models with an upper limit of a TeV under the assumed naturalness 
assumptions. The lightness of the gluino mass leads to the dominance of the gg over all other 
sparticle production processes in the GNLSP models. 



Eq.®), which has a g -»• x°g of 65 % with SUSY-HIT and 89% with ISA JET. The gluino decay 
branching ratios are discussed in further detail later. 

In the left panel of Fig. (OH) we give a comparison of the average mass of the squarks vs the 
gluino mass in GNLSP and compare it to that for the mSUGRA model. A similar analysis for 
the lightest squark mass is given on the right panel. In both cases one finds that aside from a 
very small region, the spectra from these two models do not overlap. Quite remarkably for the 
assumed naturalness conditions as described in Sec. (jHJ) the gluino mass in the GNLSP models 
has an upper limit of about a TeV, while in mSUGRA this limit extends far beyond. On the 
other hand, the upper limit on the squark masses is quite large extending to several TeV in 
both cases. It is the relative lightness of the gluino mass in the GNLSP case that enhances the 
production of the gluinos relative to all other SUSY production processes such as gq and titi 
at the LHC. This further explains the result of Fig. ([3]) and of Table (jHJ) which show that the gg 
production cross section dominates over all others. Thus in effect, in the GNLSP case, the LHC 
will become a gluino factory. Nonetheless, the production cross section for the squarks is still 
significant and their production could be detectable with efficient cuts. For example, the stop 
is usually the lightest squark in GNLSP models, so its production could be significant, and it 
can decay via t\ — > x + + b — > W + + b + x°- Since the W + has decays into / + + u, the stop 
will have some leptonic signatures while the gluino is a pure jet signal. An analysis of these 
leptonic signals requires a further dedicated analysis which is left for a future analysis. We also 
note that superficially the overlap of the few GNLSP points with the mSUGRA point might be 
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construed as the existence of a degeneracy for these parameter points. However, we need to keep 
in mind that what is plotted is an average squark mass. Further, as already discussed there is a 
significant splitting between the squark masses and the slepton masses for the mSUGRA model 
for the first two generations while there is very little splitting in this case for the GNLSP model. 
Thus at the very least the slepton-squark splittings lift any perceived degeneracy indicated in 
Fig.(j4]). We should also point out that although the production cross section of squarks and 
gluinos can be comparable for GNLSP models and for mSUGRA models, their LHC signatures 
tend to be significantly different due to the mass spectra of sparticles being rather different and 
specifically this is so since the gluino is the NLSP in the class of models we consider. 

Early LHC Discovery Prospects at Low Luminosity: Each GNLSP model is subject to the 
experimental constraints as discussed in Section([3]). We investigate the LHC signatures of 1070 
such GNLSP models with PYTHIA coupled to PGS4 [52]. Branching fractions are computed 
with SUSY-HIT and fed directly into the PYTHIA decay table via the SUSY Les Houches Accord 
(SLHA) [51] interface. For the GNLSP models, this is quite important as one must take into 
account the radiative decay g — > x°9 which can be substantial for this class of models. The 
LHC detector simulation proceeds with PGS4 with the Level 1 (LI) triggers designed to mimic 
the Compact Muon Solenoid detector (CMS) specifications [53] with the LHC detector card. 
Specifically the LI trigger level cuts that are imposed are as follows[52]: (1) inclusive isolated 
lepton (ji/e) (30 GeV); (2) lepton plus jet (20 GeV, 100 GeV); (3) isolated di-leptons (15 GeV); 
(4) di-leptons plus jet (10 GeV, 100 GeV); (5) isolated di-leptons (10 GeV); (6) isolated lepton 
plus isolated r (15 GeV, 45 GeV); (7) isolated di-tau (60 GeV); (8) inclusive isolated photon 
(80 GeV); (9) isolated di-photon (25 GeV); (10) inclusive tf T (90 GeV); (11) inclusive single-jet 
(400 GeV); (12) jet plus f T (180 GeV, 80 GeV); (13) acoplanar jet and f T with (1 < < 2) 
(100 GeV, 80 GeV) (jet, f T )\ (14) acoplanar di-jets with (A<j> < 2) (200 GeV). Muon isolation is 
controlled by employing the cleaning script in PGS4. SM backgrounds have been generated with 
QCD multi-jet production due to light quark flavors, heavy flavor jets (bb, tt), Drell-Yan, single 
Z/W production in association with quarks and gluons (Z/W+ jets), and ZZ, WZ, WW pair 
production. The standard criteria for the discovery limit are imposed, namely, that the SUSY 
signal is taken as discoverable if the number of SUSY events exceeds 5 V Asm or 10 whichever is 
larger, i.e., Nsusy > Max {5a/ Asm, 10} . We implement several classes of post trigger level cuts 
to analyze the SM background and our event samples. 

Post Trigger Level Cuts: As already noted the dominant sparticle production process at the 
LHC will be gg — > gg, and the dominant decays of the gluino are qqx° and depending on the 
particular part of the parameter space the decay g — ► gx°[5H ES, E] can also be large, and in 
fact can dominate over the 3 body decay, depending on a confluence of the following: (a) the 
lightness of the gluino, (b) the largeness of the squark masses, and (c) the neutralino mixing 
matrix. While a heavy gluino can decay into quarks and squarks and into weak gauge bosons 
W ± , Z[56], these decays are either kinematically forbidden, such as into quarks and squarks, or 
are highly suppressed for the GNLSP models we are discussing here. The above implies that 
GNLSP models will lead to a preponderance of jet signatures. In our analysis we impose three 
different sets of post trigger level cuts to optimize the signal and reduce the background from the 
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standard model processes. We classify these as: (i) Class 1 (CI), (ii) Class 2 (C2), and (iii) Class 
3 (C3) post trigger level cuts. We discuss these in some detail below and discuss the signals that 
are best detected with these three classes of cuts. Before proceeding further, we note that the 
missing transverse momentum is an important cut both as a trigger level as well as a post trigger 
level cut allowing one to increase the signal relative to the background. Thus SUSY models with 
a LSP which is massive tend to produce events at the hadron colliders with a larger missing 
energy. In order to suppress the Standard Model background, usually a large missing transverse 
momentum cut is employed. We note further that since b and b are produced in the processes 
99 ~^ QQyQ9y99 as we U as with a certain fraction in gg — > gg with subsequent decays of g (this 
latter case being the most relevant one discussed here), one has a significant number of b jets 
produced in these events. Thus b -tagging is a useful instrument in their identification. These 
features will be seen in the post trigger level cuts we discuss below. 

(i) Class 1 Post Trigger Level Cuts (CI) 

1. Electrons, and muons with Pt > 10 GeV (where Pt is the transverse momentum) and 
1 77 1 < 2.4 (where 77 is pseudorapidity) are selected. 

2. Jets with P T > 60 GeV and 1 77 j <C 3 are selected. 

3. Events with P T > 200 GeV are selected. 

Since the GNLSP models have a gluino which lies close to the LSP consistent with the relic density 
constraints, a large portion of events generated by the SUSY processes have a less energetic jet 
signal and a relatively smaller missing energy than may be observed, for example, in the case 
of stau coannihilation (see [35J for analysis of such a signature). Taking the above into account 
we modified the imposed CI post trigger level cuts in C2 and C3 which are to be discussed 
below. As already emphasized, the SUSY production cross section are dominated by gluino pair 
production and decays of the gluino. Therefore, we only select events that have at least two jets 
that pass our jet selection condition. We also investigate the effects of putting a softer missing 
Pt cut and jet Pt cut. 

(ii) Class 2 Post Trigger Level Cuts (C2) 

1. Electrons, and muons with Pt > 10 GeV and |?7| < 2.4 are selected. 

2. Jets with Pt > 50 GeV and \rj\ < 3 are selected. 

3. Events with Pt > 150 GeV are selected. 

4. Events with at least 2 jets are selected. 

The dominant SM background for GNLSP models are from QCD, Z/W+ jets, bb, and ti and 
we focus on these in Fig. ([5]). In the left panel of Fig. ([5]) we give an analysis of these backgrounds 
at the LHC with events/bin/fb^ 1 as a function of the azimuthal angle A0(jet 1; jet 2 ) between the 
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Figure 5: (Color online) Left panel: An analysis of events/bin/fb _1 as a function of the az- 
imuthal angle A0(jet 1 jet 2 ) between the two hardest jets in the GNLSP model relative to the 
SM backgrounds with C2 cuts. Right panel: H T = Ylijets distributions of GNLSP models and 
SM backgrounds with C2 cuts. These distributions act as a guide for implementing the C3 cuts 
as discussed in the text. 



two leading jets. We also show the distributions for two GNLSP model points with masses of 
400 and 500 GeV respectively. A similar analysis is presented but as a function of the total jet 
Pt (labeled Ht) in the right panel of Fig.Q, where again we also exhibit the distributions for 
two GNLSP model points. The analysis provides a good starting point for charting out the post 
trigger level cuts needed to reduce the background and enhance the signal. Thus the analysis 
suggests that one may cut out the events that have a large A0(jet x , jet 2 ) as this cut will suppress 
the QCD background due to light quark flavors, bb and ti. A veto on isolated electrons or muons 
was applied to reject the background events containing W or Z leptonic decays. A simple count- 
ing of events after applying C2 cuts reveals that the GNLSP models are nearly lepton free, as 
is the case for the bb and di-jets background. However, the standard model backgrounds due to 
Drell-Yan process and (Z/W+ jets) have (31% ~ 36%) of the total events that contain electrons 
or muons, while for ti background it is about 45%. The ZZ, WZ, WW pair production resulting 
in multi- leptonic backgrounds have an even larger percentage of leptonic events. Therefore, the 
e//i veto significantly enhances the GNLSP signals over the standard model background. This 
leads us to investigate cut class C3. 
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Figure 6: (Color online) The analysis above is with post trigger level cuts C3 and with an 
integrated luminosity of 10 fb _1 . Left panel: average Pt vs the gluino mass; Right panel: the 
discovery reach in SUSY events vs the gluino mass. 



(Hi) Class 3 Post Trigger Level Cuts (C3) 

1. Apply cut set C2. 

2. Electron or muon veto is imposed. 

3. H T = J2 jets P T >A00 GeV. 

4. The azimuthal angle A0(jet l5 jet 2 ) between jetl (the hardest jet) and jet2 (the second 
hardest jet) is chosen so that A0(jet 1; jet 2 ) < 37r/4. 

5. The azimuthal angle A0(jet l7 Pt) between jetl (the hardest jet) and Pt is chosen so that 
A^jeti,#r) >tt/2. 

6. The azimuthal angle A0(jet 2 ,-^r) between jet2 (the second hardest jet) and Pt is chosen 
so that A0(jet 2 ,^ T ) > 7r/4. 

As is indicated from the preceding discussion, jets with and without tagged b jets are important 
signals for the discovery of the GNLSP models. Another important signature is {Pt), which is 
the average magnitude of the missing transverse momentum, where the average extends over all 
events passing the cuts. We discuss now several signatures for the GNLSP models as given in 
Fig.©. The analysis of Fig.© is given at y/s = 14 TeV under post trigger level cuts C3 with an 
integrated luminosity of 10 fb^ 1 . The left panel of Fig.© exhibits the average Pt as a function 
of the gluino mass. Here one finds that essentially all of the parameter points of the GNLSP 



19 



300 400 500 600 700 800 900 
Gluino NLSP Mass (GeV) 



300 400 500 600 700 800 900 
Gluino NLSP Mass (GeV) 



Figure 7: (Color online) All Cuts and constraints as in Fig. ([6]): Left panel: the discovery reach 
with 10/fb for SUSY events with 1 tagged b-jet vs the gluino mass, and Right panel: A com- 
parison of post trigger level cuts CI and C3 for the discovery of a 4-jet signal with an integrated 
luminosity of 10fb _1 at the LHC. It is seen that the specialized cuts C3 enhance the 4J reach 
by roughly 100 GeV relative to the post trigger level cuts CI. The C3 cuts reduce the SM 
background while allowing a large n-jet signal for n > 2 with the largest signal to background 
enhancement appearing for n = 4. 



models have an average fa which is larger, and often significantly larger, than for the standard 
model case for which the average /V is found to be ~ 257 GeV under the same set of cuts. 
The right panel illustrates the ratio iVgusY / V N$m vs the gluino mass where iVsusY is the total 
number of SUSY events and Nsm is the standard model background, again with the imposition 
of post trigger level cuts C3. Here one finds that the ratio A^usy / V Asm lies above 5er discovery 
limits for gluino masses up to 800 GeV making this ratio an important channel for the discovery 
of the GNLSP models. 

The left panel of Fig.® demonstrates that the events containing a single tagged b-jet leads 
to a discoverable signal over a wide range of gluino masses. Specifically this panel gives the ratio 



^susy/ v ^sm* vs the gluino mass where Ag^gy is the number of SUSY events with 1 tagged 

b-jet and A^f* is the number of standard model events with 1 tagged b-jet under the imposition 
of class C3 cuts. This channel provides a 5a discovery for gluino masses up to about 600 GeV. 

We briefly comment on the ability to tag b-jets in this case. Naively one would expect that 
for cases in which the dominant decays are g — > with sizeable branchings in b quarks, that 
the b-jet events would be a gold plated signal. However, here the b's come out rather soft since 
the mass splitting of the (g — x°) is rather small, typically around 50 GeV (the phenomenon 
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responsible for the satisfaction of the relic density is precisely this mass split). While the ability 
to tag the b's is possible, it is indeed more difficult than the canonical situation seen on the 
Hyperbolic Branch of REWSB [ID] (for recent work on b-tagging analyses see [121 [57J [37]). 

Nonetheless, as discussed above the signal for singly tagged b-jets is strong enough that it 
can be a useful one. The right panel of Fig. ([TJ) we compare the 4-jet discovery limits under 
CI and C3 cuts. Here it is found that the SM background drops by an appreciable amount as 
one goes from CI cuts to C3 cuts, allowing one to extend the discovery limit in 4-jet channel 
by over 100 GeV with C3 cuts relative to imposing the CI cuts in this channel. Amongst the 
classes of n-jet events, we find this channel to be the most enhanced when passing from CI to 
C3. Combining all the channels analyzed in Fig.( l6f7I) one finds that gluino masses up to 800 GeV 
are discoverable by this technique for a GNLSP model with just 10 fb _1 of integrated luminosity. 
Thus the validity of the GNLSP models can be tested with first data from the LHC. 



8 Dark matter detection in the GNLSP models 



Many of the GNLSP model points have an LSP neutralino which is bino like. However, there 
is also a significant set of models that have an LSP with large higgsino components (many of 
which are at high tan/? but this is not exclusively so). The large tan/? parameter points are 
easily spotted by examining the tan/3 vs Ao/mo plot in Fig.®. The GNLSP models have 
important implications for the direct detection of neutralino dark matter. An analysis of the 
spin independent and spin dependent cross sections in dark matter experiments is given in 
Fig.® implemented with MicrOMEGAs. Included are published limits from the ZEPLIN-III 
experiment [58] . the first five-tower CDMS data[59] and the XENON 10 results [60]. Projected 
limits (indicated by *) from CDMS and LUX are also shown [6Tj . The a si vs m^o analysis 
shows some interesting results. First one finds that there are a class of GNLSP models which 
are beginning to be constrained by the direct detection dark matter experiments. These models 
would also produce large gg production cross sections and will be easily visible at the LHC. 
However, interestingly there is another class of low mass neutralino (and low mass gluino) models 
which have rather small spin independent cross sections which are outside the reach of the direct 
detection experiments in the foreseeable future. Nonetheless these models would still lead to 
rather large gg production cross sections and hence will be visible at the LHC. Thus the LHC 
can detect many of the GNLSP models which are most likely inaccessible to the direct detection 
dark matter experiments. On the flip side, for the large set of low mass GNLSP models one still 
has the possibility of a very light neutralino (gluino) with a sizeable higgsino component and 
consequently a large spin independent cross section. Thus if a light gluino is indeed indicated 
early on at the LHC, it may also provide a hint of the size of the dark matter signal in direct 
detection of dark matter. We note in passing that a plot of asi vs m g looks very similar to 
the left panel of Fig.® as the gluino and neutralino mass are related at the low scale via 
trig = (1 + A^o)m^o. 
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Figure 8: (Color online) Left panel: An exhibition of the spin independent cross section asi as a 
function of the neutralino mass. It is seen that there are a large number of models corresponding 
to asi i n the range 10 -46 cm 2 and below, which would be inaccessible in direct dark matter 
searches in the foreseeable future. However, many of these models especially those with low 
values of m^o (and hence of rrig) would be discoverable at the LHC even with low luminosity. 
Right panel (1): The explicit scaling relation between m^o and m g for the GNLSP models (also 
shown purely for visual reference is a line representing m^o = rrig) . Right panel (2): An analysis 
of asi vs the light Higgs boson mass illustrating that a large portion of NUSP13 (light blue) has 
a Higgs boson near 120 GeV while NUSP14 is given in dark (magenta). 



We note that in the GNLSP class of models, the direct annihilations of x°X° into electron 
positron pairs is helicity suppressed. For example, we obtain (av)^o^o^ e + e - ~ 5 x 10~ 30 cm 3 /s 
(at v/c = .002) for the model of Tabled while the self annihilation into W + W~ are equally 
small. Annihilations into rr are found to be the largest (this particular model has {av)^o^Q_ Tf ~ 
2 x 10~ 27 cm 3 /s). Thus a significant boost will be needed to explain the recent cosmic ray 
excess [52]. 



9 GNLSP Benchmarks 



It is useful to give benchmarks for the three GNLSP models A, B, and C discussed in Sec.(T2]). 
As mentioned above, we go beyond the perturbative calculation, and include the Sommerfeld 
enhancement of the cross section in these benchmarks. In Table we give benchmarks for model 
GNLSPa- The benchmarks are chosen to exhibit a significant diversity in the input values. 
Some of the low lying spectrum as well as the spin independent cross section asi{x°p) an d 
the spin dependent cross section <Jsd{x°p) corresponding to Tabled are exhibited in Table([7|). 
The analysis of Table (j7j) shows a variation over two orders of magnitude for asi{x°p)- As 
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discussed already the variation arises due to changes in the higgsino vs the bino component of 
the LSP. Similar benchmarks for model GNLSPb are given in Table (jHJ) and the corresponding 
light sparticle spectrum and the corresponding spin independent cross section <Jsi{x°p) an d the 
spin dependent cross section <Jsd{x°p) are exhibited in Table([9]). Finally, a similar analysis for 
the model GNLSP C is given in Table (PBJ and Table dHJ). 

10 Conclusions 

In the above we have given an analysis of a class of models with nonuniversalities which lead 
to a gluino as the NLSP. Several important observations emerge from this analysis which have 
bearing on the observation of sparticles at the LHC Perhaps the most important of these is 
that if the gluino is the NLSP, then the gg — > gg cross section at the LHC dominates over all 
others in the GNLSP models. The dominance of the gluino production and the fact that the gg 
production cross sections are large implies that the observation of supersymmetry via the gluino 
production can occur with the first data from the LHC. It is found that the dominant signal 
of the gluino NLSP model are multi-jets, tagged b-jets and missing energy and it is possible 
to devise post trigger level cuts which discriminate these models above the standard model 
backgrounds. Such cuts which reduce the background and enhance the signal to the background 
ratio were devised and implemented in this paper. We note also that the inverse of the LHC 
process, namely gg — > gg, is largely responsible for the satisfaction of the relic density when 
the neutralino and the gluino coannihilate. Further, the analysis of GNLSP models reveals that 
there exists a significant region of the parameter space in these models where the neutralino has 
a large higgsino content, and the neutralino-proton spin independent scattering cross section is 
sizeable and can be probed with the current experimental sensitivity, and sensitivities that would 
be achievable in future experiments. However, there are also other regions of the parameter space 
where the neutralino is mostly bino like and in this case the spin independent cross sections can 
fall well below the current experimental sensitivity, and well below the sensitivity that would 
be achievable in the near future experiments. Interestingly, the bino cases, though difficult to 
discover in dark matter experiments, can be accessible at the LHC since the gluino mass in 
these models lies within the reach of the LHC even at low luminosities. Another aspect of 
the GNLSP models was also discussed which relates to the compressed nature of the sfermion 
mass spectrum relative to the case of universal gaugino masses. Here the sleptons in the first 
two generations could be almost degenerate in mass and often even heavier than their squark 
counterparts. Several benchmarks for the GNLSP models were also given to facilitate further 
work. It was also pointed out that a test of the GNLSP models can be done with just 10 fb _1 of 
integrated luminosity and thus it is one of the models that can be checked with the early data at 
the LHC. Finally, we emphasize once again that if the gluino is an NLSP then the production of 
gluinos will dominate all other sparticle production making the LHC effectively a gluino factory. 
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23 



Appendix A: Gaugino mass sum rules for GUTS with nonuni- 
versalities 



Each of the mass ratios listed in Table© gives two gaugino mass relations at the GUT scale. 
Of these one is 'unstable' to the inclusion of a singlet F term breaking while the other one is 
'stable' and remains valid when one includes a singlet F term breaking along with the non- 
singlet breaking. Below we list only the 'stable' mass relations. The mass relations are labeled 
numerically (l)-(23) as in Table©. They ar^| 

(1) + (7) + (10) : -5Mi + 3M 2 + 2M 3 
(3) : -Mi + 9M 2 - 8M 3 
(5) : -5M X + 9M 2 - 4M 3 
(13) : Mi + M 2 - 2M 3 
(16) : -5Mi + M 2 + 4M 3 
(20) : 5Mi + 7M 2 - 12M 3 

In addition to the above, models (6), (8), (11), (12), (15), (17), (22) in Table© satisfy the 
relation M 2 = M 3 , while model (19) in Table© satisfies the relation Mi = M 3 , and model 
(21) in Table© satisfies the relation Mi = M 2 . These mass relations would be appropriately 
modified at low scales by the renormalization group evolution. Thus at the one loop level the 
mass relations at the electroweak scale are 

a 3 (0) 



0. 


(2): 


Mi + 3M 2 - 4M 3 = 0, 


0. 


(4): 


5Mi + 3M 2 - 8M 3 = 0, 


0. 


(9): 


5Mi + M 2 - 6M 3 = 0, 


0. 


(14): 


-25M X + 9M 2 + 16M 3 = 


0. 


(18): 


5Mi + 11M 2 - 16M 3 = 0, 


0. 


(23): 


-5Mi + 9M 2 - 4M 3 = 0. 



M V(Q) + ^#A.f 2 (Q) 



-M 3 (Q) = 0, 



(16) 



OL\{Q) OC 2 (Q) Oi 3 (Q)' 

where r is the ratio as given by Table©. Assuming one can determine with accuracy the gaugino 
masses, these mass relations can be a useful indicator of the specific F type breaking and hence 
of the type of nonuniversal SUGRA model one has at the GUT scale. 



Appendix B: An analytic analysis of quasi-degeneracy of 
LSP and GNLSP and of the sfermion mass compression 

Here we give a one loop analysis of how the quasi degeneracy of the LSP and of the GNLSP comes 
about and then discuss a similar phenomenon for the squarks and the sleptons for the first two 
generations. We begin with the gaugino masses which at the GUT scale obey the nonuniversality 
condition of Eq.© which we restate in more compact notation as follows: For model i we have 

M[ i] = 7i l) m 1/2 , M 2 W = 7 «mi /2 , M« = 7 <V 1/2) (17) 

9 The stability of the gaugino mass relations to the inclusion of singlet breaking is easily seen by noting that 
the sum of the co-efficients of Mi, and M3 in each of the mass relations in Eq. (fl6)) vanishes. 



24 



where 



7 ? } = (1 + a iai ), 7? = (1 + bioa), 7 f = (1 + (18) 

With the above boundary conditions, the gaugino masses at the electroweak scale t = ln(M G /Q 2 ), 
where M G is the GUT scale and Q is the electroweak scale, are given by 

M®{t)=>y®a a {t)m x , 2 , a = 1,2, 3. (19) 

Here a=l,2,3 correspond to the U(l), SU{2) and SU(3) gauge groups, and a a (t) are the corre- 
sponding fine structure constants for these groups at the electroweak scale. For the mSUGRA 
case, 7a ^ = 1, and one simply has that at one loop 

Mi(t) : M 2 (t) : M 3 (t) = «i(t) : a 2 (t) : a 3 {t). (20) 

Using the experimental values of the gauge coupling constants at the electroweak scale one finds 
the three gaugino masses roughly in the ratio 1 :~ 2 :~ 5 — 6. In this case the gaugino masses 
are split in a very significant way. However, in the presence of nonuniversalities the ratios will 
be modified in a very different way. Thus with the inclusion of the modifications j a , the three 
gaugino masses at the electroweak scale will be roughly in the ratio 71 :~ 272 :~ (5 — 6)73. It 
is clear then that the choice 73/71 ~ (1/5 — 1/6) will make masses of gaugino 1 and of gaugino 
3 almost degenerate. Of course, more realistically there would be mixings between the gauginos 
and the higgsinos and the mass eigenstates will be admixtures of these. Thus the mass relation of 
the lightest neutralino and of the gluino will be more complicated. Still the above approximation 
may roughly hold when the neutralino is mostly a Bino. 



Nonuniversalities also enter in the masses for the squarks and for the sleptons. For the first 
two generation down squarks one finds 

m LW = m o + m 2 di + a G [-f 3 + -f 2 + 3q/iK/ 2 + (-3 + 3 sin2 B w) M l cos(2/?), 

m\ R {t) = ml + ml + a G [^f 3 + ^fi]m 2 1/2 - i sin 2 9 W M 2 Z cos(2/5), (21) 

where 

/. = *., /.w = s( 1 -(r+W) 1 (22) 

where (3 a = b a a a (0), & a (0) = a a /4ii, and b a = (33/5,1,-3) for the gauge groups U(l), SU{2) 
and SU(3). For the case of the first two generations of charged leptons one has 

m\ L {t) = ml + m 2 ei + a G [-f 2 + ^fi]m 2 1/2 + (-- + sin 2 9 W )M 2 cos(2/3), 

m l R (t) = m l + m l + ^OL G hm\ j2 - sin 2 6 W M 2 Z cos(2/5). (23) 
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Using the above we find that 

\ [ m LW + - (™L(*) + m lM] = m l - <i + «g[^7 3 2 / 3 - ^7i/i]- (24) 

It is now easily seen that the last brace in Eq. (1241) can vanish or even turn negative by appropriate 
choice of 73 vs 71. The above situation leads to a near equality of the squark and of the slepton 
masses discussed in the text of the paper. The /i parameter also has a strong dependence on 
nonuniversalities. This is seen by examining the relation that determines /1 2 , i.e., 

V? = (m 2 Hl ~ m 2 H2 tan (3 2 ) (tan (3 2 - l)" 1 - -M 2 Z + Afi 2 , (25) 
where A/1 2 is the loop correction. Thus m^, and m# 2 are sensitive to nonuniversalities since 

m 2 Hi =m 2 + a G Q/ 2 (t) + rn 2 1/2 , (26) 

where the nonuniversalities enter via the / functions. Similarly m# 2 is given by 



m 



2 _ _2 ~ 
h 2 



m{ /2 e(t) + A m m 1/2 f(t) + m 2 (h(t) - k(t)A 2 ), (27) 

where the tilde functions e and / are modified due to nonuniversality while the functions h(t) 
and k(t) are not unaffected (for definitions of these functions see the first paper of [15] and 
the references therein which also gives a more detailed discussion of this topic). The depen- 
dence on nonuniversalties is more complicated in this case because of the coupling with the top 
quark. Specifically, one can derive the following relation which gives the explicit dependence on 
nonuniversalities 

|^ = (t 2 - l)- l {m\ /2 g' a - t 2 {m\ /2 e' a + A m m 1/2 f a )) + (28) 

where g' a = J^-, g — 5g(|/2 + e ' a = an d /„ = q- One can make a semi-quantitative 

estimate of the dependence of /i 2 on nonuniversalities from above. We note, however, that \i 
does not enter sensitively in the sum rule for the first two generation of squarks and sleptons and 
thus an estimate of the compression of the sfermion spectrum in the first two generations can be 
made without estimate of the \x parameter. 



Appendix C: Sample GNLSP benchmarks 



We provide here sample model points for the GNLSP class of models. Each model point obeys 
experimental constraints as discussed in the text. We use here SuSpect 2.41 coupled to Mi- 
crOMEGAs (MO) 2.2.CPC along with an independant code which agrees with MO at the per- 
tuabtive level but accounts for the non-perturbative effects discussed in the text. Similar model 
points may be obtained with other spectrum calculators coupled to MO [see, for example, Ta- 
blefll]) where one such comparison is given, which, however, is only at the perturbative level]. 
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GNSLP 


m (GeV) 


mi/2 (GeV) 


A (GeV) 


tan /3 






GNSLPai 


2949 


692 


3658 


35 


0.566 


-0.847 


GNSLPa2 


2706 


783 


4408 


37 


0.560 


-0.839 


GNSLPas 


2529 


946 


3873 


41 


0.560 


-0.837 


GNSLP A4 


2967 


910 


5114 


27 


0.557 


-0.834 


GNSLPas 


2574 


1058 


4197 


42 


0.557 


-0.833 


GNSLP A6 


2821 


1019 


5050 


20 


0.554 


-0.830 


GNSLP A7 


3008 


1252 


-3241 


27 


0.558 


-0.836 


GNSLPas 


2746 


1265 


5186 


17 


0.551 


-0.824 



Table 6: GNLSPa benchmarks: These models produce the correct relic density with SuSpect 
2.41 coupled to MO 2.2.CPC including the non-perturbative corrections . Benchmarks here have 
(tth 2 ) x o G (0.100,0.130). Here m t (pole)/ GeV =170.9 throughout. 



GNSLP 

Model 


m h 
(GeV) 


(GeV) 


m-g 
(GeV) 


(GeV) 


m t 
(GeV) 


m A \h\ T 2 ° 

(GeV)j- 


°si(x°p) 
(Pb) 


o"sd(x°p) 
(Pb) 


GNSLPai 
GNSLP A2 


117 
117 


285 
337 


343 
387 


343 
869 


1560 
1225 


2130 | 0.241 
1952 0.004 


3.0xl0~ 8 
2.7xl0~ 10 


8.0xl0" 5 
2.8xl0~ 7 


GNSLPas 
GNSLP A4 


116 
117 


398 
399 


456 
454 


480 
1080 


1190 
1276 


1540 | 0.135 
2724 0.003 


2.3xl0~ 8 
1.3xl0~ 10 


2.4xl0~ 5 
l.lxl0~ 7 


GNSLPas 
GNSLP A6 


116 
117 


447 
448 


510 
507 


531 
1064 


1161 
1149 


1514 | 0.131 
2886 0.003 


2.2xl0~ 8 
2.1xl0~ 10 


1.9xl0~ 5 
1.4X10" 7 


GNSLP A7 
GNSLPas 


120 
116 


551 
557 


618 
624 


645 
970 


1332 
1032 


2718 | 0.108 
2960 0.007 


1.3xl0~ 8 
6.7xl0~ 10 


LOxlO" 5 
3.3xl0- 7 



Table 7: GNLSPa properties: The variation in the spin independent cross section over two 
orders of magnitude arises due to variations in the higgsino vs the bino component of the LSP. 
The higgsino fraction has been defined by hj 1 ^ — |t| 2 + |<^| 2 ) where the normalized LSP mass 
eigenstate is x° = ab + (3w + 7/11 + 5h 2 as in the notation defined in Sec.(|4]). 
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GNSLP 


m (GeV) 


mi/a (GeV) 


A (GeV) 


tan (5 


6 2 


5s 


GNSLPbi 


2421 


736 


3414 


51 


0.000 


-0.841 


GNSLP B2 


3406 


734 


4655 


35 


0.000 


-0.848 


GNSLPbs 


2890 


945 


-2977 


47 


0.000 


-0.844 


GNSLP B4 


3772 


988 


5894 


46 


0.000 


-0.837 


GNSLP B5 


2857 


1158 


-2631 


24 


0.000 


-0.842 


GNSLP B6 


2943 


1142 


5006 


12 


0.000 


-0.831 


GNSLP B7 


3188 


1376 


-2479 


7 


0.000 


-0.837 


GNSLPbs 


2659 


1380 


5028 


37 


0.000 


-0.825 



Table 8: GNLSPb benchmarks: As in TableQ the displayed models produce the correct 
relic density. Benchmarks here have (Qh 2 ) x o e (0.100,0.120). 



GNSLP 

Model 


m h 
(GeV) 


(GeV) 


m- g 
(GeV) 


(GeV) 


(GeV) 


I !, Frac 

mA\rii 2 
(GeV)|- 


°si(x°p) 
(pb) 


°sd(x°p) 
(Pb) 


GNSLPbi 
GNSLP B2 


116 
119 


313 
316 


364 
364 


557 
576 


1223 
1734 


409 | 0.009 
2520 | 0.003 


2.4xl0" 8 
1.4xl0- 10 


1.2xl0~ 6 
2.0xl0~ 7 


GNSLPbs 
GNSLP B4 


118 
119 


417 
428 


484 
488 


582 
779 


1507 
1818 


665 | 0.035 
1756 | 0.002 


2.0xl0~ 8 
1.0xl0~ 10 


5.2xl0~ 6 
5.6xl0~ 8 


GNSLP B5 
GNSLP B6 


119 
117 


502 
502 


566 
565 


583 
948 


1475 
1299 


2610 | 0.140 
3177 0.002 


1.8xl0~ 8 
1.7xl0- 10 


1.7xl0~ 5 
6.8xl0~ 8 


GNSLP B7 
GNSLPbs 


117 
116 


601 

598 


672 
669 


701 
1051 


1676 
1121 


3382 | 0.113 
2031 0.004 


1.8xl0~ 8 
4.3xl0^ 10 


7.4xl0~ 6 
1.6xl0~ 7 



Table 9: GNLSPb properties: The table gives an analysis similar to that of Table ((Zj) for 
GNLSPb models. 
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GNSLP 


m (GeV) 


mi/2 (GeV) 


A (GeV) 


tan p 


s 2 


5 3 


GNSLPci 


1604 


450 


2035 


49 


-0.317 


-0.852 


GNSLP C2 


2119 


696 


2860 


44 


0.291 


-0.845 


GNSLP C3 


2443 


948 


2823 


11 


-0.227 


-0.842 


GNSLP C4 


3850 


1111 


4388 


9 


0.209 


-0.840 


GNSLP C5 


2599 


1270 


3656 


14 


-0.009 


-0.833 


GNSLP C6 


2458 


1479 


5045 


35 


0.462 


-0.823 


GMSLPc7 


2087 


453 


2359 


21 


0.292 


-0.862 


GNSLP C8 


1958 


674 


2950 


22 


0.299 


-0.843 


GNSLPcg 


3874 


1098 


4455 


8 


-0.018 


-0.839 


GNSLP C10 


2543 


1337 


4188 


49 


-0.147 


-0.826 


GNSLP C11 


3288 


1431 


5879 


15 


0.703 


-0.822 


GNSLP C12 


3942 


1755 


5995 


36 


0.080 


-0.825 



Table 10: GNLSPc benchmarks: A sample of benchmarks in GNLSPc in a random distri- 
bution in 62 and £3. Many of the models listed above have a substantial higgsino component and 
part of the parameter space would be accessible to current and future experiments for the direct 
detection of dark matter. Benchmarks here have (Qh 2 ) x o £ (0.100,0.138). 



GNSLP 

Model 



(GeV) 



m x° 
(GeV) 



m-g 
(GeV) 



(GeV) 



m t 
(GeV) 



VIA I h\*2 ° 

(GeV)j- 



<7Sl(X°P) 

(Pb) 



3.6xl0- 8 
1.4xl0~ 8 
1.3xl0~ 8 
1.4xl0~ 8 
2.9xl0^ 8 
2.0xlQ- 8 



°sd(x p) 
(pb) 



9.0xl0~ 6 
1.6xl0~ 5 
1.2xl0~ 5 
l.OxlO" 5 
1.9xl0" 5 
9.7xl0" 6 



GNSLPci 
GNSLP C2 
GNSLP C3 
GNSLP C4 
GNSLPcs 
GNSLPce 



113 
115 
115 
117 
116 
115 



185 
291 
404 
481 
542 
632 



227 
338 
461 
546 
610 
704 



228 
423 
508 
587 
609 
716 



830 

1073 

1315 

2078 

1320 

830 



516 

1065 

2488 

3939 

2677 

2042 



0.022 
0.056 
0.075 
0.098 
0.200 
0.130 



5.1xl0" 9 
1.2xl0~ 9 
l.lxl0~ 9 
5.4xl0~ 9 
1.7xl0~ 9 
3.7xl0^ 9 



2.2xl0" 5 
1.8xl0~ 6 
6.0xl0~ 7 
1.3xl0~ 6 
8.0xl0~ 7 
1.6xl0- 6 



GNSLPct 

GNSLPcs 

GNSLP C9 

GNSLPcio 

GNSLPcn 

GNSLP C i2 



115 
114 
117 
116 
117 
119 



189 
289 
483 
576 
634 
769 



223 
333 
546 
648 
708 
849 



333 
581 
816 
797 
909 
948 



1114 
924 
2088 
1242 
1310 
1965 



1938 
1887 
3996 
956 
3558 
2909 



0.047 
0.012 
0.010 
0.018 
0.015 
0.031 



Table 11: GNLSPc properties: A display of a part of the sparticle mass spectrum consisting 
of light Higgs and CP odd Higgs masses, and masses of the LSP, NLSP, light chargino, and 
light stop along with spin independent and dependent cross sections. Models shown here include 
those with higgsino like LSPs and as well as those with mixed higgsino and bino LSPs, and LSPs 
which are mostly bino. The horizontal line emphasizes separation of the osij pb from O(10~ 8 ) 
to O(10~ 9 ) . 
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